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Section 1 
SUMMARY 


Proton-induced pulse height distributions (PHDs) in Si:XX detectors were studied analytically 
and experimentally. In addition, a preliminary design for a flight experiment to characterize 
the response of Si: XX detectors to the trapped proton environment and verify PHD models 
was developed. 

PHDs were computed for two orbit altitudes (200 and 900 km) for a variety of shielding 
configurations (1 to 19.3 gm/cm 2 ). Most of the proton-induced pulses have amplitudes less 
than about 3.5 x 10 5 e-h pairs. Shielding has a small effect on the shape of the PHDs (for 
practical cases). The primary effect of shielding is to reduce the total number of pulses 
produced. 

Proton-induced PHDs in a Si:Sb focal plane array bombarded by a unidirectional 67-MeV 
beam were measured. The maximum pulse height recorded was 6 x 10 5 pairs. The 
distribution had two peaks: the larger peak corresponded to 3.8 x 10 5 pairs and the smaller 
peak to 1 .2 x 10 5 pairs. The maximum pulse height and the larger peak are within a factor of 
two of predicted values. The low-energy peak was not expected, but is believed to be an 
artifact of inefficient charge collection in the detector. 

The planned flight experiment will be conducted on a shuttle flight. The Lockheed Missiles & 
Space Company, Inc., helium extended life dewar (HELD) will be used to provide the 
required cryogenic environment for the detector. Two bulk Si:Sb arrays and two Si:As 
impurity band conduction arrays will be tested. The tests will be conducted while the shuttle 
passes through the South Atlantic Anomaly. PHDs will be recorded and responsivity changes 
tracked. This experiment will provide a new database on proton-induced PHDs, compare two 
infrared-detector technologies in a space environment, and provide the data necessary to 
validate PHD modeling. 
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Section 2 
INTRODUCTION 

The major tasks in this program were theoretical and experimental investigations of the 
effects of trapped protons on the operation of infrared detectors in the South Atlantic 
Anomaly (SAA) and the design of a flight experiment to enlarge the database on proton 
effects on infrared focal plane arrays (FPAs) and validate proton pulse height models. The 
theoretical efforts were development of a PHD code and application of the NOVICE code to 
the computation of PHDs. In the experimental portion of the program, Si:Sb detectors were 
irradiated with 67-MeV protons, and radiation-induced signals were collected and analyzed. 
Before discussing the results of these studies, we present a brief phenomenological 
discussion of the problem. The flight experiment design will be discussed in a later section. 

Ionizing particle radiation, such as that encountered in the trapped radiation belts and in 
strategic engagement, produces random pulses in the outputs of optical and infrared 
detectors. These pulses can be a significant source of noise. The number of noise electrons 
N generated in the detector in time r is 

N = Wc 

where n is the average pulse generation rate and \ Is the number of hole-electron pairs 
produced per interaction in the detector, n is assumed to be Poisson distributed and A. is a 
random variable. The probability of a pulse of a given magnitude being generated (the PHD) 
is a function of the particle type and energy, the detector geometry, and the type and quantity 
of material surrounding the detector. The pulse generation rate (event rate) depends on 
these quantities and on the particle flux density. 

The pulse event rate due to charged particles, e.g., electrons, protons, and alphas, is a 
straightforward though difficult calculation. After accounting for the particles’ energy losses in 
the detector surroundings, the PHD for an isotropic flux is obtained by a convolution of the 


2-1 


X90823_FJ_2 


LMSC-F279280 


detector chord length distribution with the particle flux energy loss distribution. Some useful 
order-of-magnitude estimate, however, can be obtained from simple phenomenological 
arguments. 

The proton-induced event rate, Er, in a single detector element is 

E, = 0 p (/w + ti + wt)/2 events/s 

where <f> p is the isotropic proton flux density and I, w, and t are the detector dimensions. Each 
event produces a pulse of electron hole pairs N p 

N p = (C/e)dE/dX (pairs) 

where C is a chord length, e is the energy required to create an electron hole pair (3.65 eV in 
Si), and dE/dX is the energy loss rate of the proton. The average pulse height can be 
estimated by inserting average values for the quantities in the above equation. Typical 
dimensions of a bulk extrinsic Si detector pixel are 75 x 75 x 300 p 3 . The peak proton flux 
density in a 300-km circular orbit (28° inclination) is about 200/cm 2 -s. This flux density will 
produce an event rate of 0.05 events/s per pixel. A significant fraction of the total proton flux 
falls in the 50- to 200-MeV range. The average proton energy loss rate in that energy interval 
is 1 .2 x 10 3 eV/pm. The average chord length of the detector is about 75 pm. The average 
value of the PHD should be about 2.5 x 1 0 4 pairs. Noise levels at the detector array output of 
less than 200 electrons (rms) can be achieved. Thus, during a 500-s integration time, one 
expects about 25 events which exceed the detector noise level by about 1000 times. These 
events, due to their frequency and magnitude, significantly reduce the performance of the 
detector. 
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Section 3 
RESULTS 

In this section we present the results of both our analytical work on calculating pulse height 
distributions and the results of the cyclotron tests on an Si:Sb FPA. 

3.1 PULSE HEIGHT DISTRIBUTION CALCULATIONS 

The calculation began with the computation of the chord length distribution of a detector. The 
starting point was a paper by Bradford ( J. Appl. Phvs. . Vol. 50, 1979, p. 3799), which gives 
a general expression for the sum distribution of chord lengths in a rectangular parallelepiped. 
The chord length distribution is the negative derivative of the sum distribution with respect to 
the chord length. Several mistakes were found in Bradford's paper and were corrected (see 
Appendix A). The final results for the sum distribution may be determined by analytical 
evaluation of integrals or by numerical integration of intermediate expressions. Both 
techniques were used to reproduce one of Bradford’s results. The numerical integration 
technique gave the better results. The failure of the analytical technique is believed to be 
related to roundoff errors since the expressions are very long and complicated. 

At this point in the program we acquired the NOVICE code, a mature Monte Carlo transport 
code which had been used elsewhere to compute pulse height distributions in detectors. 
Further development of our own pulse height distribution code was halted for reasons of 
economy. 

The NOVICE computer program is basically a particle transport and shielding calculation 
code. Our particular use of the many features of the NOVICE program has been to calculate 
pulse height distributions (PHDs) for Si detectors within a spherical shield exposed to an 
isotropic proton flux. Proton flux densities appropriate to the SAA were used. The energy 
distribution of the proton flux was separately calculated using the AP8MIC model. NOVICE 
performs a Monte Carlo simulation of sensor irradiation with randomly selected incident 
proton energy and angle of approach to the detector. 
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Through a series of input data sections the user defines the Si sensor and shield 
configurations. The detector investigated was a centrally located detector in a 58 x 62 Si 
detector array. The array was 300-pm thick and each pixel was 75 x 75 p 2 . The detector 
array was In bump bonded to a Si readout chip (500-p.m thick). The whole array was 
mounted on a Cu heat sink. The array and heat sink were surrounded by a spherical shell 
1-cm thick; the effect of shield material was investigated. The user defines the desired output 
of the program: in this case, the range of the calculated pulse height spectrum. Pulse heights 
are given in energy deposited per event. To convert energy (in eV) to electron hole pairs, 
divide the pulse height by 3.65. 

Pulse height data are accumulated and output to a special file at the end of the calculations. 
This file is available for display or further analysis. Output data are in units of the number of 
events (energy spikes) within each designated subrange of the pulse height spectrum, 
normalized to the incident radiation flux. To obtain the count density, we divide each count 
output of the NOVICE program by the width of its respective energy interval. 

The uncertainty of the results due to taking a finite number of incident particle samples in the 
Monte Carlo simulation process are displayed in the PHD figures by plotting the mean value 
and the +/- one standard deviation value. Higher energy values are subject to greater 
uncertainty than lower energy values for a given scenario, because fewer of the randomly 
selected incident particles yield pulses in the higher energy range. 

We first investigated the effect of proton energy and shielding density on PHDs from some 
simple spectra. The results of calculations are easily related to phenomenological models 
and provide a good check on the operation of NOVICE. Three quasi-monochromatic uniform 
isotropic proton energy spectra were used: 48 to 52 MeV, 68 to 72 MeV and 98 to 102 MeV. 
A 1-cm shield of density either 1 gm/cm 3 or 2.7 gm/cm 3 surrounded the focal plane array. 
The resultant PHDs are shown in Fig. 3-1. The lower energy input spectra produce larger 
pulses than the highest energy spectrum. The maximum energy a 98-MeV proton can 
deposit in the detector is about 0.5 to 1 MeV (taking into account the energy loss in the 
1-cm-thick shield). Since the proton energy loss rate decreases with increasing energy, this 
shift in PHD is expected. The lower portion of the PHD is not as strongly dependent on proton 
energy: however, the steep portion of the distribution does move toward larger pulse 
amplitudes with decreasing proton energy. 
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Fig. 3-1 Normalized PHD, Density 1 Shield. Narrow-Band Incident Spectra 

Increasing the density of the shielding material should have the same effect as lowering the 
proton energy. This is observed (see Figs. 3-2, 3-3, and 3-4). The peak energies of the 
energy-degraded spectra are 39, 45, and 82 MeV. We note that the energy-degraded 68- to 
72-MeV spectrum produces a PHD very similar to the undegraded 48- to 52-MeV spectrum 
(compare Figs. 3-1 and 3-3). 

The mean chord length c of a rectangular parallelepiped is 

2twh 

C Sw +£h + wh 

The mean chord length of a 75- x 75- x 300-p 3 detector is 67 \i. The energies deposited by 
the mean energy proton in each of the three spectra are: 1.5 x 10 5 , 1.2 x 10 5 , and 8.9 x 
10 4 eV for the 48- to 52-, 68- to 72-, and 98- to 102-MeV spectra, respectively. These mean 
energies and their trends with energy compare well with the result of the NOVICE calculations 
(Table 3-1). The higher pulse count for the p = 2.7 shield relative to the p = 1 shield is 
believed to be due to a softening of the proton energy spectrum by the denser shield. 
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Fig. 3-2 Normalized PHDs for 48- to 52-MeV Input Spectrum 



Fig. 3-3 Normalized PHDs for 68- to 72-MeV Incident Spectrum 
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Fig. 3-4 Normalized PHD, 98- to 102-MeV Incident Spectrum, Shield Densities 1.0, 2.7 


Table 3-1 SUMMARY OF PHD CHARACTERISTICS FOR QUASI- 
MONOCHROMATIC PROTON FLUXES 


Energy Band 
(MeV) 

Mean of PHD 
(MeV) 

Shield Density 

Variance of PHD 
Shield Density 

Relative Number 
of Pulses* 

P - 1 

B 

P = 1 

■SB 

48-52 

0.19 

0.72 

0.028 

0.37 

1.4 

68-72 

0.14 

0.21 

0.014 

0.033 

1.6 

98-102 

0.10 

0.12 

0.005 

0.009 

1.2 


•Relative Number of Pulses = 


Total Number of Pulses for p = 2.7 Shield 
Total Number of Pulses for p * 1 Shield 
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PHDs for two trapped proton spectra were calculated for different shield densities. The PHDs 
from these broader energy spectra are more complex and difficult to interpret. The proton 
spectra for 300-km and 900-km circular orbits with 28® inclinations were computed using 
AP8MIC. These two energy spectra are summarized in Table 3-2. 


PHDs were computed for four shield configurations: 1-cm shield (p = 1, 2.7, 19.3) and an 
AI-Ta-AI combination (2.5 x 10' 2 cm Al + 0.1 cm Ta + 0.23 cm Al). 1 This latter shield 
configuration was chosen because it effectively removes trapped electrons. The minimum 
shielding configuration (1 cm of 1 gm/cm 3 material) is well below minimum shielding 
configuration for any practical case. The 1 cm of 2.7 gm/cm 3 material corresponds to the 
expected shielding configuration in the planned flight experiment. The effects of slight 
perturbation in the proton energy spectra were also investigated by uniformly shifting the 


Table 3-2 ORBITAL PROTON ENERGY DISTRIBUTIONS 


Proton Energy 
(MeV) 

Relative Number Protons/MeV 
in Each Band 


300 km 

900 km 

4-6 

1.13 x 10- 2 

8.4 x 10' 3 

6-10 

1.15 x 10' 2 

7.4 x 10‘ 3 

10 - 15 

8.2 x 10- 3 

7.7 x 10' 3 

15 - 30 

6.23 x 10- 3 

6.33 x 10- 3 

30 - 50 

5.65 x 10- 3 

5.00 x 10‘ 3 

50 - 70 

5.57 x 10- 3 

4.88 x 10- 3 

70 - 90 

4.50 x 10‘ 3 

4.22 x 10- 3 

90 - 110 

4.47 x 10‘ 3 

3.93 x 10- 3 

110 - 150 

3.68 x 10‘ 3 

3.31 x 10' 3 

150 - 200 

2.15 x 10- 3 

2.26 x 10' 3 

200 - 250 

1.32 x 10~ 3 

1.52 x 10- 3 

250 - 300 

6.96 x 10- 4 

9.78 x 10- 4 

300 - 350 

3.76 x 10- 4 

6.3 x 10- 4 

350 - 400 

2.06 x 10' 4 

4.08 x 10" 4 

400 - 450 

1.86 x IQ’ 4 

7.54 x IQ’ 4 
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input spectra by ± 0.25 MeV. Characteristics of the computed PHDs are summarized in 
Table 3-3. 

The general trends of all the PHDs computed are well illustrated by Figs. 3-5. 3-6. 3-7, and 
3-8 /or the 900-km input spectrum and the four shielding configurations. The shielding effects 
at low-pulse amplitudes can more easily be seen in Fig. 3-9. 

Increasing shield density primarily affects the number of pulses produced. According to 
Table 3-3, the total number of pulses (per unit incident flux density) decreases by a factor of 
1.5 as the shield density increases from 1 to 19.3 gm/cm 3 . This reduction occurs over the 
entire PHD. 

The broad peak around 5 MeV is observed in all the PHDs computed using broadband 
energy input spectra. This peak is not the result of a distortion in the proton energy spectrum 
introduced by the shielding computation in NOVICE. 

This was verified by computing the PHD for the case of no shielding and a truncated energy 
distribution. The high end of the energy distribution was truncated by energy corresponding 
to the energy absorption in the shield. This PHD displayed all the features of the PHD shown. 

A series of diagnostic PHDs was calculated using different energy input spectra and no 
shielding. The spectra were uniform over the following energy bands: 2 to 10, 2 to 30, 2 to 60, 
and 2 to 200 MeV. The evolution of the PHDs with increasing proton energy can be seen in 
Fig. 3-10. Consider first the 2- to 10-MeV spectrum case. The highest energy pulse is 
produced by a proton whose range is just equal to the longest chord length of the detector 
(320 mm) . Higher-energy protons have a lower-energy ion rate and therefore deposit less 
energy along this chord. 

Lower-energy pulses will be produced by both higher-energy protons and lower-energy 
protons. For example, a 10-MeV proton (whose range is about 770 p.m) will deposit 2.6 MeV 
along the longest chord. A 2-MeV proton will deposit all its energy in 25 p.m. Thus, the low 
end of the PHD can be attributed to both low- and high-energy protons. 

Figure 3-1 1 shows the chord length sum distribution and chord probability density for a 75- x 
75- x 300- p 3 detector. Note that 50% of the chords have lengths greater than 75 |im (the 
range of a 2.5-MeV proton). Thus, low-energy protons will be effectively absorbed, producing 
pulses in the MeV range. 
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Fig. 3-9 Normalized PHD, 300-km Spectrum, 1-cm Shield 



Fig. 3-10 Normalized PHD, No Shield. Various Flat Spectra 




X80823_FJ_3 


LMSC-F279280 



CHORD LENGTH (pim) 


Fig. 3-11 Chord Length Distribution, 75- x 75- x 300-p.m Volume 

As higher-energy protons are added to the spectrum, the low-amplitude pulse density 
increases. These protons have ranges far in excess of any chord length and therefore 
produce small pulses. The PHD for the 2- to 200-MeV spectrum is very similar to the PHDs in 
Figs. 3-5 through 3-8, including a shallow trough at about 1.25 MeV. The depth of the trough 
in Fig. 3-8 may be due to the detailed shape of the proton distribution after passing through 
the shielding material. Recall that there was no shielding in the cases shown in Fig. 3-10. 

In summary, it is clear from these calculations that shielding reduces the total number of 
pulses produced in the detector. Large amounts of shielding (e.g., 19.3 gm/cm 2 ) do not 
significantly reduce the number of pulses. This is due to the fact that most of the pulses are 
the result of interactions with high-energy protons which cannot be shielded against. 

Shielding and fluctuations in the trapped proton energy distribution can affect detector noise. 
Mandel (Brit. J. AddI. Phvs .. Vol. 10, 1959, p. 222) developed a formalism for calculating the 
noise associated with image fluctuations. The formalism is applicable to proton-induced 
noise processes In detectors. Consider n incident particles, each of which creates X. 
secondaries (X. is a random variable) . We associate n with the incident proton and X. with the 
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number of electron hole pairs in a proton-generated pulse in the detector. The variance A p 
in the number of electron-hole pairs produced is 

A V = X 2 A^n 4- nA 2 * 
where the variances all have the form 


A 2 /* = n 2 -n 2 etc. 

This holds for any distribution of n and \. Adding shielding changes all the quantities on the 
RHS. In addition, a perturbation in the trapped photon energy spectrum will affect all the 
parameters in the equation. A fluctuating energy spectrum will add a new noise term. We can 
use the data in Table 3-3 to investigate the effect of shielding and energy distribution 
fluctuation on noise. 

In making these comparisons, it is convenient to equate the total number of pulses to n and to 
divide A*jl by n 0 . the total number of pulses for the p = 1 shielding case. In this way a direct 
comparison of each can be made. Since n is Poisson distributed, (An) 2 = n, the above 
equation becomes 


no n 0 L J 

The results for the 900-km orbit are summarized in Table 3-4. 

For a constant input energy distribution, shielding always reduces the noise. A small shift in 
the energy spectrum toward lower energies can increase the noise by as much as a factor of 
two. The noise is then very sensitive to the actual proton distribution, and fluctuations in the 
spectrum will produce a new noise term. 

3.2 FPA IRRADIATION TESTS 

Proton irradiations of a Si:Sb FPA were made at the University of California at Davis (UCD) 
cyclotron facility. A photograph of the 58 x 62 array is shown in Fig. 3-12. The detector 
array is In-bump bonded to a switched field-effect transistor (FET) multiplexer readout circuit. 
A schematic diagram of the circuit is shown in Fig. 3-13. The FPA was contained in the Ames 
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Table 3-4 IMPACT OF SHIELDING AND ENERGY SPECTRUM ON NOISE 


Shielding 

Shifted Energy 
Distributions 

P : 

= 1 

P = 

2.7 

P = 

i 

H 

19.3 

Al Ta Al 

+0.25 

-0.25 

+0.25 

, -0.25 

n/no 

i 

0.83 

0.71 

0.88 

0.91 

0.91 

0.79 

1.2 

A 7 ^ 

0.038 

0.025 

0.024 

0.029 

0.028 

0.05 

0.021 

0.049 

no 










“Cryotran” test dewar (Fig. 3-14). A block diagram of the test electronics is shown in 
Fig. 3-15. Pre-experiment testing and evaluation of the FPA were performed by the NASA 
Ames personnel. Both the dark can and the vacuum jacket had thin Al windows. A 1 -in. -thick 
Al plate with a 0.75-in. hole was placed immediately in front of the dewar vacuum window to 
define the proton beam. The proton beam traveled about 10 in. in air before entering the 
dewar. Two plastic scintillation counters were positioned in the air space. By using 
coincidence counting, an accurate determination of the proton flux density could be made. 
The proton beam energy was 67 MeV, a value chosen for reasons of beam control and 
stability. 

One limitation of ground-based experimentation is that only monoenergetic or nearly 
monoenergetic particle beams are available. The pulse height distribution from such a beam 
is not, in itself, interesting. As an example of such an “uninteresting” distribution, we show 
the PHD generated by a uniformly distributed 40- to 60-MeV proton beam (Fig. 3-16). The 
PHD should consist of a single well-defined peak corresponding to the energy deposited 
along the length (interelectrode spacing) of the detector (for our normal incidence irradia- 
tions). Since the pulse peak is well defined, deviations from this ideal are readily observed. 
These deviations may be due to proton interactions in the FPA surrounding which produce 
secondary sources of ionizing radiation or charge collection anomalies in the detector array. 
Examples of secondary ionizing radiation fluxes are: K and L-shell x rays, bremsstrahlung, 
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BUCK AND WHITE PHOTOGRAPH 



Fig. 3-12 58 x 62 Si:Sb Infrared Focal Plane Array 
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Fig. 3-13 Schematic Diagram of Multiplexer Unit Cell 


VACUUM SHROUD ATTACHED 
TO CRYOSTAT VIA FOUR 
TAPPED-THROUGH 6-32 


SCREWS ON 8 11/1 6-in. -DIA. VACUUM 

BOLT CIRCLE SHROUD 
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Fig. 3-16 Pulse Height Relative Distribution (Unidirectional) 

and y photons produced in the detector surroundings. Charge-collection anomalies include 
incomplete charge collection and charge sharing (i.e., blooming). 

The objectives of these experiments were: 

(1) Determination of the PHD for a well-defined (energy and angle of incidence) 
proton beam 

(2) Observation of blooming due to the high density of deposited charge in a 
pixel 

Coincidence pulses due to a proton traveling through two detectors would not be observed in 
our experiments. 

The FPA used was fabricated by the Santa Barbara Research Center. The detectors were 
300-jx thick with 75- x 75-p. 2 pixel areas. The Cryotran dewar, FPA, and electronics were 
supplied by NASA Ames. 

The first cyclotron run produced no useful PHD data. The first proton hit caused the array to 
hang up at the supply voltage. With the beam off, the array could be reset. The first proton hit 
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after reset would initiate hang-up again. Post-irradiation diagnostics revealed that a combina- 
tion of a low voltage on V g0 and a faulty FPA ground provided the necessary conditions for a 
proton-induced pulse to initiate hang-up. 

In the second test, the output of a pixel well away from the edges of the array was monitored. 
A centrally located pixel was chosen to eliminate complications due to nonuniform collection 
fields which may be present at the edges of the array. In addition, it requires much less time 
to read out a single pixel than it does to read out the whole array. The objective was to collect 
as many data in as short a time as possible. From other studies we have found that in order to 
obtain statistically valid PHD, about 10 4 pulses need to be counted. Typically, 100 pulse 
height channels would be used. If 10 4 counts were uniformly distributed over all the channels, 
then a rms error of 10% in each channel is expected. In the case of a nonuniformly distributed 
PHD, the peak of the distribution will be defined with greater precision than the wings of the 
distribution. 

Proton flux densities of 1700 and 4000 p/cm 2 -s were used. Detection integration times were 
16 or 32 ms. Table 3-5 is a matrix of the test conditions. 

The number of events expected E p is 


Ep — (pp Aq Tj Nj 

where 

<j> p = proton flux density 
A d = detector area 

Table 3-5 TEST CONDITIONS 


Proton Flux 

Integration Time 

Samples 
Recorded at 

Density 

16 ms 

32 ms 

16 ms 

32 ms 

1700 

— 

V 

0 

2.8 x 10 4 

4000 

V 0 


2 x 10 4 

4.9 x 10 4 
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t i = integration time 

N| = number of integration samples collected 

Table 3-6 summarizes the expected number of events expected for each combination of flux 
density and integration time. 

The total number of events expected was about 510; the observed number of events was 
~ 5300, over 10 times the expected value. The resulting pulse height distribution is shown in 
Fig. 3-17. The distribution is quasi-exponential with a peak pulse height of about 0.27 V. In 
this run, the amplifier gain, G, was 16 and the sense capacitance, C, was 7 x 10‘ 14 f. The 
conversion of pulse voltage to pairs is 

CV 

N p = — P alrs 

where e is the magnitude of the electronic charge. 

The peak pulse contained 7.4 x 10 3 pairs and corresponds to an energy deposited 
of 2.7 x 10 4 eV (3.65 N p ). The largest pulse at 1.75 V contained 4.6 x 10 4 pairs 
(1.7 x 10 5 eV deposited). The energy loss rate for a 67-MeV proton in Si (p = 2.33 g/cm 3 ) 
is 1.9 x 10 7 eV/cm (“Studies in Penetration of Charged Particles in Matter,” Barkas and 
Berger, National Acad. Sci., Nuclear Science Series, Report No. 39, 1964). The expected 
energy deposited in a 300-|im-thick detector is about 5.6 x 10 5 eV. The entire PHD is well 
below this expected value. The observed PHD is anomalous with respect to both number of 
events and the magnitude of the events. The excess number of pulses could be due to 


Table 3-6 EXPECTED NUMBER OF 
PROTON EVENTS 



r, 

16 ms 

32 ms 

1700 

86 

72 

4000 


350 
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Fig. 3-17 Proton PHD (E p = 67 MeV) 

(1) Inaccurate dosimetry 

(2) Extraneous sources of ionizing radiation produced by protons in the detector 
surroundings 

(3) Collection of charge from adjacent pixels 

(4) Events in the multiplexer 

Because of the type of dosimetry used, we assume at this time that (1) is not a viable 
explanation. 

Proton interactions in the array surroundings can result in the emission of ionizing radiation 
which may produce pulses in the detectors in addition to the pulses due to direct proton hits. 
These interactions can produce 

• Bremsstrahlung 

• X rays from K- and L-shell ionization events 

• Gamma photons from nuclear interactions 
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The bremsstrahlung production rate is inversely proportional to the mass of the incident 
particle. Because of this dependence, the bremsstrahlung production rate by protons is 
always negligible compared to the production rate by electrons (R. B. Evans, The Atpmip 
Nucleus . McGraw-Hill Inc., 1955, p. 600). The impact of trapped electron bremsstrahlung on 
the Shuttle Infrared Telescope Facility (SIRTF) focal plane array was investigated by Autio 
(NASA Report, Dec. 1975). He concluded that electron bremsstrahlung-lnduced pulse event 
rates were negligible compared to proton event rates in the detector. We concluded that 
proton bremsstrahlung-included events would not be observed. 

Protons can ionize K- and L-shell electrons, resulting in the emission of x rays from the 
detector surroundings. A single proton can produce a shower of x rays which would affect a 
number of detectors in the array. Cross-section data for x-ray production by high-energy 
protons are not available. However, since the cross sections peak in the 1- to 2-MeV range 
and decrease at higher energies (E. Merzbacher and H. W. Lewis, Encyclopedia of Physics , 

Springer-Verlag, Berlin, 1958, p. 167), use of lower energy data will result in a worst case 

-20 

estimate of the effect. The maximum cross section for x-ray production in Al is 3 x 10 
cm 2 . Since the absorption coefficient of K- and L-shell x rays is on the order of 10' 2 cm, only 
those protons interacting in the detectors' immediate surroundings will generate x rays which 
can reach the detector. The number of x-ray photons, I, per proton is 

I = Nct< 5 photons/proton 

where N is the atomic density of the material surrounding the detectors (Al) , o is the cross 
section for x-ray production in Al, and 8 is the x-ray absorption length in Al. We find that 
I = 30. In our tests, the detector array was about 2 cm behind an Al cap. The maximum x-ray 
fluence at the array is about 1 .2 x-ray photons/cm 2 -proton. The Al K(alpha) emission occurs 
at 1.55 keV. Each x-ray photon will generate about 400 electron hole pairs. This effect could 
result in an increased number of small (relative to a direct proton hit) pulses in a detector 
output. 

Cross sections for (p, gamma) reactions were calculated at UCD. The gamma photons 
produced in these interactions have energies in the 1- to 5- MeV range. These photons have 
meanfree paths of about 4 cm. Thus, gamma photons produced in the collimator can reach 
the detector array. The collimator will be the largest source of gamma photons since it totally 
absorbs a portion of the beam. To estimate the worst-case gamma flux density at the 


3-22 


X90823_FJ_3 


LMSC-F279280 


detector, we assume that the proton beam is totally absorbed in the collimator and that no 
photons are absorbed by the dewar walls. The (p, gamma) cross section in Al at 67 MeV is 
4.4 x 10" 27 cm 2 . The gamma flux density (G) at the array is 

G = Ncr t 0 P A/( 27 rR 2 ) gamma photons/cm 2 /s 

where N is the atomic density of Al, o is the cross section, t is the thickness of the collimator 
(2.5 cm), 4> P is the proton flux density, A is the beam area (8 cm 2 ), and R is the collimator 
array separation (18 cm). Using these values, we find that G = 2,6 x 10" 6 . Clearly, this 
effect is negligible. 

Option (3) is very attractive. The centrally located pixel is surrounded by eight pixels. We 
expect about 510 events in each pixel; therefore, in all nine pixels we would expect a total of 
about 4600 counts, close to the observed number. Pulses from adjacent pixels would be 
smaller than those originating in the central pixel, due to recombination. If we take the 
analogy of optical cross-talk, then these pulses should be about only 1% as large as the 
pulses from the central pixel. Thus, we would expect a PHD that is skewed toward small 
pulse heights. 

The peak in the observed PHD corresponds to the energy deposited by a 67-MeV proton in 
14 pm of Si. This distance is not unlike the distance of the sum of a p-n junction depletion 
width and a minority carrier diffusion length in a Si readout circuit. Thus, proton events in the 
readout circuit could be cited as the source of pulses observed (begging the question of what 
became of pulses due to events in the detector). However, since the nodal areas in the 
readout circuit are an order of magnitude smaller than detector areas, we would expect that 
this contribution to the PHD would be smaller than the contribution due to events in the 
detectors. This effect is then inconsistent with the observed event rate. It should also be 
noted that in Co 60 irradiations of FPAs, few or none of the observed pulses can be attributed 
to interactions in the readout circuit. 

In the third cyclotron irradiation test, the PHD was constructed from events in all pixels. 
Approximately 100 frames were recorded. These frames contained about 10 4 proton events. 
In these runs, detector bias was -5 V. The free field pixel signals were less than 200 mV. 
The data (thresholded at 2 V) are presented in Figs. 3-18 and 3-19. The data in the first 
figure were taken from the first 50 frames; the data in the second figure from the last 50 
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frames. Both PHDs terminate at about 11 V and have peaks at about 7 V. The gain in all 
cases is 8 and the sense capacitance is 7 x 10 -14 F. The 7 V peak corresponds to the 
generation of 3.8 x 10 5 pairs (1.4-MeV deposited). The maximum value pulse contains 
about 6 x 10 5 pairs (2.2-MeV deposited). The PHDs have a low voltage peak at 2.25 V 
(1.2 x 10 5 pairs). Attempts to generate a PHD from a unidirectional monochromatic proton 
beam have not been successful. Only two points were computed after in excess 16 K Monte 
Carlo histories. The "triangular” PHD (Fig. 3-16) comes closest to representing the 
conditions existing when the data of Figs. 3-18 and 3-19 were taken. The peak in Fig. 3-16 is 
at an energy of 0.73 MeV (2 x 10 5 pairs); the pulse height at a frequency 100 times lower 
than the peak frequency is 1 .4 MeV (3.8 x 10 5 pairs). These values are about a factor of two 
lower than those observed. However, since the computation did not exactly simulate the 
experimental conditions, the argument is considered satisfactory. 

The simulation predicts a very symmetric PHD, in contrast to the observed spectra, which 
exhibit an additional peak at low energies. Other investigations using monochromatic 
gamma-ray sources have observed PHDs which are skewed toward low pulse heights 
(private communication, Mark McKelvey, NASA Ames). Figure 3-20 shows the response of a 



BIN CENTER VOLTAGE 

Fig. 3-18 PHD From Unidirectional 67-MeV Proton Beam (Sum of Frames 13-69) 


3-24 





X90823_FJ_3 


LMSC-F279280 


Si:Sb detector array to 60 keV, 241 Am photons. The maximum voltage pulse at 0.575 V is 
due to a 1 .6 x 10 4 pair pulse. This corresponds to 60 keV, the maximum energy the gamma 
ray can deposit. The preponderance of low-amplitude pulses suggest that the charge 
collection efficiency is incomplete and that most pixels have low charge-collection efficien- 
cies. The same phenomenon may be active in the proton experiments producing the low 
energy peak. 

Initial analysis of the data Indicates that blooming does occur. However, blooming does not 
appear to extend beyond nearest neighbor pixels. In summary, the observed proton PHDs 
are in general agreement with computer predictions. The low energy peak is likely due to 
charge-collection inefficiencies in the detector rather than to an anomalous proton interac- 
tion. 
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Section 4 

FLIGHT EXPERIMENT CONCEPT 

In this section we describe a flight experiment whose objectives are to establish a database 
on proton-induced PHDs in Si: XX detectors, compare the response of bulk and impurity band 
conduction (IBC) detectors to proton bombardment, and validate proton-induced PHD 
modeling in Si:XX focal plane arrays. This flight experiment study was prepared within a 
framework based on the needs, goals, and limitations of the Outreach Program. 

4.1 JUSTIFICATION FOR FLIGHT EXPERIMENT 

SIRTF and some DoD missions plan to use Si:XX and other IR detectors. Depending on the 
orbit chosen, these detectors may pass through the trapped radiation belts. The trapped 
electron flux can be excluded with modest shielding: the proton flux cannot be shielded 
against. These protons will generate charge pulses on the order of 10 5 pairs in the detectors 
and will significantly degrade the signal-to-noise ratio of the detection system. In order to 
develop algorithms to reduce the effects of these pulses, their pulse height distribution must 
be known. Codes have been developed to compute the PHDs. However, the results of the 
codes have never been tested. Ground testing cannot accurately simulate the omnidirec- 
tional, energy-distributed proton flux encountered in space. A flight experiment provides the 
opportunity to test these codes. The validated PHD model then can be used with confidence 
to predict PHD in other detectors, e.g., InSb and HgCdTe. 

Bulk Si: XX detectors have been in use for a number of years. IBC detectors are now available 
in limited numbers and are expected to replace bulk Si:XX detectors. This flight experiment 
would also provide an opportunity to compare the response of detectors of the two 
technologies in identical radiation environments. This experiment will enlarge our database 
on proton radiation effects in Si:XX detectors, provide needed PHD data, validate PHD 
modeling, and make a direct comparison of the responses of bulk and IBC detectors. 

4.2 SUMMARY 

The experiment will be performed on a shuttle flight. A 5-day experiment is planned. The 
experiment will require minimal attention by the flight crew (power on and power off only). 
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The experiment will be contained in the Lockheed Missiles & Space Company, Inc. (LMSC) 
project Helium Extended Life Dewar (HELD). The dewar will be mounted on one of the new 
extension capabilities of the Space Transport System (STS), the HITCHHIKER-M (HH-M). 
Some of the electronics will be mounted on the dewar. The remaining electronics will be 
stored in a sealed HH-G canister. Data will be collected during and after passage through the 
SAA and then down-linked using the medium-rate KU-band downlink. Data analysis will be 
performed on the ground by LMSC personnel. Another LMSC experiment (the Microsphere 
Investigation) will share the instrument space of the HELD. Since the focal plane experiment 
requires only 3.5 x 3.5 in., there is the possibility that a third experiment could be housed in 
the HELD instrument chamber. 

Before describing the flight experiment hardware, we present a brief outline of the experiment 
protocol. The objectives of this experiment are to establish a database on trapped proton 
PHDs in two types of extrinsic Si detectors, validate pulse height prediction modeling, and 
assess long-term degradation effects. To obtain the largest possible database, the PHDs will 
be derived from the outputs of all the detectors on the FPA. 

After achieving a stable orbit, the experiment electronics will be powered up. The experiment 
will contain its own computer to control measurement sequence. First, clear environment 
responsivity and noise measurements will be made (clear environment meaning outside the 
SAA. The FPAs will be operated continuously during passage through the SAA and all frames 
will be stored. After the third pass through the SAA, another clear environment check of 
responsivity and noise will be made. This agenda will be followed for the remainder of the 
experiment. Data will be downlinked as necessary. Periodically, after the post-SAA respon- 
sivity checks, the FPAs will be cycled to 15 K or higher to anneal out any total dose effects. 

4.3 FLIGHT EXPERIMENT REQUIREMENTS 
4.3.1 Data Requirements 

The prime objective of this experiment is to gather statistically meaningful data on trapped 
proton-induced PHDs in infrared detectors of interest to NASA programs. The primary need is 
data fidelity. We must record a large number of pulses in order to achieve a statistically 
meaningful sample. Data and focal plane requirements appropriate to Si:XX detectors will be 
used in the following discussion. Similar requirements will apply if other detector technolo- 
gies, e.g.. InSb or Ge:Ga, are incorporated in the experiment. 
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A 300 km orbit was chosen for this discussion. The results will not be significantly different if 
the experiment were performed in a 900 km orbit. The AP8MIC model was used to calculate 
proton flux density and energy distribution. The total number of protons of all energies 
encountered in the SAA in 1 day is 2.38 x 10 5 P/cm 2 . Typical effective areas of bulk Si:XX 
detectors are about 2.5 x 10' 4 cm 2 (75-M.m x 75-pm x 300-jim detector). A single pixel 
would experience about 60 proton hits/day. An IBC detector with the same optical area would 
receive about four times fewer proton hits. Since our experiment will last only a few days, this 
approach cannot provide the 10 4 data samples we require. If. however, we read all the FPA 
pixels, then we will obtain about 2 x 10 5 data points/day in a bulk technology FPA and about 
4 x 10 4 data points/day in an IBC technology array (each array will contain 58 x 62 
elements). Thus, in order to achieve the needed database, the PHD will be constructed from 
the data from all the pixels on the FPA. This is not an optimum experimental design since 
edge pixels may have a different response to protons than do interior pixels (due to different 
fringing electric fields around the pixels) . However, these external pixels constitute less than 
10% of the total pixel population. We believe that this tradeoff between number of data points 
and absolute fidelity is required for a viable experiment. 

Total dose effects due to protron bombardment are not of major interest in this experiment. 
These experiments can be adequately performed on the ground. The total dose in the course 
of this mission will be small. A single 50-MeV proton deposits about 3 x 10” 3 r/pixel hit. 
Therefore, a single pixel will accumulate about 0.5 r per day. If the flight lasts 5 days, the total 
dose will be about 2.5 r. Total dose effects at 1 r have been observed in ground tests. Thus, 
these measuremets may be significant. 

4.3.2 FPA Characteristics and Requirements 

Extrinsic silicon infrared detectors are strong candidates for IR astronomy applications such 
as SIRTF. Such detectors can be doped with In, Bi, As, Ga, or Sb to achieve various spectral 
band coverage. In particular, Sb-doped Si detectors (Si:Sb) can cover the spectral band 
from 15 to 31 pm with a peak response at 27 p.m. When mated to a suitable readout circuit, a 
responsivity >3 A/W and a total read noise of less than 130 rms electrons has been 
demonstrated at an operating temperature of 8 K 1 . These detectors require large thick- 
nesses to achieve adequate quantum efficiency. This makes them particularly sensitive to 
radiation effects such as proton radiation. Protons passing through the devices create 
electron-hole pairs as they interact with the constituent atoms in the lattice, producing 
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unwanted noise and false signals in the devices. Also, proton interactions within the readout 
circuitry can cause fluctuations in the detector bias, which affects the responsivity of the 
detectors. These effects are clearly undesirable for IR astronomy applications. 

Currently under development are IBC extrinsic Si IR detectors that have several potential 
advantages over standard photoconductive (PC) devices. The impurity band devices are 
much thinner than the corresponding PC devices and thus are less sensitive to radiation 
effects. Also, the IBC devices do not display some of the other negative effects observed in 
standard PC devices such as the hook effect and long-time constants. The most promising 
devices tested to date are Si: As IBC devices that can effectively cover the 12- to 24-p.m 
spectral band. 

The Si:Sb detectors and the Si:As IBC devices fabricated in a 58 x 62 staring array format 
have been chosen for the flight experiment. The Si:Sb devices reach the longest cutoff 
wavelength of all the Si detectors. To date, no Si:Sb IBC devices have been demonstrated. 
The Si As detectors provide a look at the devices that will eventually replace the PC devices. 

4.3.2. 1 Description of Detectors . 

Sl;Sb . Si:Sb detectors are fabricated from Czochralski-grown Si doped to approximately 
4 x 10 15 cm' 3 with Sb. The usual concerns regarding residual shallow impurities in extrinsic 
Si detectors are not important in this case, since the Sb level is shallower than the residual 
impurities B and P. The only concern here is that the residual B concentration be as low as 
possible so that it does not compensate the Sb. Double-sided polished wafers with 
thicknesses of 300 *un are used to fabricate the detector arrays. Both surfaces are 
ion-implanted with P, followed by an Al deposition to provide ohmic contacts. The Al layers 
are delineated such that the illuminated side has windows to allow the IR radiation to pass into 
the material, while individual detector contacts are delineated on the opposite side. 

Si:As . Si: As IBC detectors are fabricated from epitaxially grown material on undoped Si 
substrates. A heavily doped layer is grown on an undoped Si substrate followed by an 
undoped epitaxial blocking layer. Detectors are fabricated in a manner similar to that used for 
the bulk PC Si:Sb devices. 

One of the main advantages of the IBC device is that the absorbing layer is only 15 to 20 ^im 
thick. Normally the heavy doping would cause impurity band current to flow in the device, 
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thus creating a shunt current that prevents this impurity band current from flowing. The result 
is a device that is less sensitive to radiation effects, has a fundamental J2 improvement in 
signal-to-noise ratio, has an extended wavelength response, and does not display the 
detrimental time constant effects of an equivalent bulk PC extrinsic Si detector. 


4.3.2.2 Description of Readout Circuit . The Si:Sb and Si:As IBC devices can be In bump 
bonded to a switched-FET readout/multiplexer that has been used by several groups in the IR 
astronomical community. 1 ' 3 A detailed description of the operation of the Si:Sb with this 
readout circuit and a characterization of its performance for space-based astronomy are 
described elsewhere. 1,4 This direct readout (DRO) is an n-channel Si MOSFET switching 
device with two parallel outputs. Unit cell pairs, which can be accessed in any order, are 
selected by means of an X,Y addressing scheme using row and column address clocks. 
Each unit cell of the DRO has a reset switch and source follower for each detector. In 
addition, there is an output source follower for each of the two outputs to drive the output 
lines. 


4 . 3 . 2.3 Focal Plane Requirements . The FPA requirements for the flight experiment are 
displayed in Table 4-1. We have chosen an existing design for the array format and 
readout/multiplexer to minimize cost and schedule impacts on the flight experiment. As 
pointed out in the above discussion, a Si:Sb starring array of this format has been 
characterized for space-based IR astronomy applications. The operating temperature and 
FPA noise requirements listed in Table 4-1 are derived from these characterization data. 
These are reasonable requirements for the Si:As IBC devices as well. The integration time 
and uniformity requirements reflect the real needs of a space-based IR astronomy applica- 
tion. Operability is specified at a much lower level than would be required for real 
applications, but for the flight experiment we do not require 100% pixel operability. Lowering 
the operability requirements will help minimize costs of procuring the FPAs. The responsivity 
is chosen so that we get reasonable signal-to-noise for this experiment, but not so high as to 
drive the cost of procuring the arrays. Temperature stability and uniformity requirements are 
derived from measured temperature dependence of noise and responsivity of these types of 
devices. For an astronomy application in which one is trying to image a faint and distant 
object, the requirements would be more stringent, probably ±0.1 K. In our experiment, we 
need only enough stability to observe proton events that create large signals. 
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Table 4-1 FPA REQUIREMENTS FOR PROTON RADIATION FLIGHT EXPERIMENT 


DETECTOR TYPES 

Si:Sb bulk, Si:As IBC 

DETECTOR FORMAT 

58 x 62 

PIXEL SPACING 

75 mm (center to center) 

OPERATING TEMPERATURE 

8 K 

RESPONSIVITY AT 2 V APPLIED BIAS 
Si:Sb 

>1 A/W 

Si:As IBC 

>1 A/W 

INTEGRATION TIME 

0.2-500 S 

FPA NOISE/INTEGRATION TIME 

100 e - rms 

OPERABILITY 

>60% 

UNIFORMITY OF OPERATING PIXELS 
Responsivity 

±10% 

Noise 

±10% 

TEMPERATURE STABILITY 

±0.5% 

TEMPERATURE UNIFORMITY OVER FPA 

±0.5 K 

NOISE LEVEL AT RECORDING INSTRUMENT 

200 e~ rms 

POWER DISSIPATION PER PIXEL 

<5 pW 

CABLE BETWEEN FPA AND DRIVE ELECTRONICS 
Noise Issues 

For lengths greater than 0.6 m, line drivers will be required 

Type of cable 

35-wire ribbon cable/array 

Low Impedance Ground 

1 /Array 

COLD FINGER AREA REQUIRED TO MOUNT 2 FPAs 

5 x 10 cm 


We propose using two of each type of FPA on this experiment. This plan will increase both 
reliability and the volume of data collected. With a slight loss in reliability, two other detector 
technologies could easily be accommodated in the HELD experiment tunnel. Candidate 
detector technologies are: InSb at 4 K and Ge:Ga operated at 2 K. 
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The clocks and biases required to run the FPA will be mounted on the exterior of the cryostat. 
FPA-to-electronics cabling will be heat sunk at an intermediate temperature of 50 K to reduce 
heat losses. In addition to the operating electronics, an infrared radiation source will be 
included in the cryostat for responsivlty measurements. A joule heated thin-film Si resistor will 
be used for the IR signal source. 

4.3.2.4 References . 

1. McKelvey, M. E,, et al., "Characterization of Direct Readout Silicon:Sb and 
SilicomGa Infrared Detector Arrays for Space-Based Astronomy." Proc. Soc. 
Photo-Opt. Instrum. Eng., Vol. 868, 1987, p. 73 

2. McLean, I. C., Proc. Workshop on Ground-Based Observations With Infrared 
Array Detectors, Hilo, Hawaii, Wynn-Williams. C. G., and Becklin, E. E., eds., 
p. 180, 1987 

3. Fowler, A. M., et al., Proc. Workshop on Ground-Based Observations With 
Infrared Array Detectors, Hilo, Hawaii, Wynn-Williams, C. G., and Becklin, E. E., 
eds., p. 197, 1987 

4. Orias, G., and Campbell, D., “Development of 30-p.m Extrinsic Silicon Multiplexed 
Infrared Detector Array," NASA Contractor Report 177, 446, 1986 

4.4 EXPERIMENT IMPLEMENTATION 
4.4.1 Mechanical/Cryogenic System 

In order to provide the required thermal environment for the FPAs, the experiment will be 
integrated into the instrument tunnel of the LMSC project HELD. HELD is an independent 
development project that has demonstrated LMSC’s superfluid helium capabilities. Some 
details of HELD are shown in Fig. 4-1. 

HELD provides a convenient design that allows easy access to experiments installed in its 
instrument well and would allow experiments that require SFHe cooling to be integrated into a 
larger experiment package that would require only one flight. The purpose of this dewar is to 
demonstrate long-term storage of SFHe in orbit for a variety of purposes. By using 
technological advances such as improved multilayer insulation and a new support strut 
system called the Passive Orbital Disconnect System (PODS), a multiyear SFHe storage 
system has been developed. 
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Fig. 4-1 HELD 

The support system for HELD is PODS, an advanced, low-conductance support system that 
has been extensively tested on the ground but not yet tested in space. Concern about PODS 
and the possible dynamic coupling between sloshing liquid helium (LHe) has led to the 
suggestion that a flight demonstration/investigation be conducted since this coupling cannot 
be measured or fully investigated on the ground. Therefore, in addition to serving as a 
testbed for the FPAs, using HELD will allow characterization of the dynamic properties of 
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PODS in a low-gravity environment. Results of the characterization of the PODS in-orbit 
dynamic response can be applied to other liquid-cryogen cooling/storage systems. HELD 
and PODS have been designed and built and are currently undergoing extensive thermal and 
structural testing. (In-space testing of HELD is described in "Integrated Cryogenic Experi- 
ment (ICE) Microsphere Investigation NASA Contractor Report 177538.) 

Many specific technologies will be demonstrated by using HELD in support of experiments 
such as 

(1) PODS 

(2) Improved multilayer insulation (MLI) 

(3) Low-conductance fill line 

(4) Porous-plug vent 

(5) PODS-supported vapor-cooled shields 

(6) Simplified prelaunch servicing 

The experiment package is designed to be mounted inside the superfluid cooled tank of 
HELD. The design of HELD permits the package to be mounted and demounted by removing 
the vapor-cooled shields and folding back some insulation. The experiment chamber will 
share a common vacuum with HELD. Instrumentation wires can be routed down PODS to 
minimize heat loads. Hermetically sealed connectors in the vacuum shell mounting ring of 
HELD provide the interfacing with flight electronics. The FPA drive electronics and signal A/D 
converters will be attached to the vacuum shell mounting ring. The remaining electronics will 
be housed in an HH-G canister. 

The HELD instrument tunnel provides a 2 K environment. The focal planes will be operated at 
8 K. In order to reduce heater power to achieve this temperature, the FPAs will be mounted 
on an Al platform thermally isolated from the 2 K sink with fiberglass/epoxy stand-offs. The 
FPA intermediate temperature stage is shown in Fig. 4-2. 

Another experiment chosen for the outreach program will share the instrument tunnel with the 
FPA experiment. This experiment will characterize the thermal properties of glass micro- 
spheres. The FPA intermediate temperature stage will be mounted on the microsphere 
experiment canister. The integration of the two experiments is shown in Fig. 4-3. 
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The four FPAs require 1 40 electrical leads (which may have resistances as high as 200 to) , 
and four leads (which must have resistances less than 20 to). Using 44-gage Ni wire for the 
high-resistance leads and 44 gage Au wire for the low-resistance leads, we compute a heat 
load of about 22 mW. The 1500 to resistance heaters will add another 8-mW heat load. The 
largest heat load will be due to focal planes themselves. The four focal planes will consume 
about 72 mW. The total thermal load of the experiment will then be less than 110 mW. The 
heat load of both experiments may be as high as 350 mW. According to Fig. 4-4, this heat 
load can be tolerated for a 5-day mission. 


One advantage of HELD is that it offers simplified prelaunch servicing by having a normal 
boiling point (NBP) helium guard tank on the first vapor-cooled shield. The purpose of this 
guard tank is to allow the SFHe tank to be filled and valved off prior to launch. The NBP tank is 
then filled and provides a 4.2 K guard shield around the SFHe tank, resulting in extremely low 
heat rates to the SFHe. The NBP tank is allowed to vent during ground hold and can 
repeatedly be refilled as needed prior to launch. 



Fig. 4-4 Effect of Experiment Heat Rate on HELD Lifetime 
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4.4.2 Experiment Electronics 

4.4.2.1 Description . The experiment electronics consists of flight- and ground-support 
hardware elements. The flight electronics provides control signals for the focal plane arrays, 
stores and encodes the experimental data, and interfaces with the orbiter's avionics 
package to receive uplink commands and transmit downlink data. The ground support 
electronics (CGSE) provides the necessary interfaces to the customer/carrier ground 
support equipment (CCGSE), including those to support the low- and medium-rate ground 
support equipment (LRGSE and MRGSE) data links. The CGSE is a workstation that enables 
the experimenter to monitor the health of the experiment’s hardware as well as sample and 
control the acquisition and telemetering of data. 

4. 4.2.2 Installation . The flight electronics consists of two assemblies, the focal plane drive 
electronics (FPDE) and the data interface unit (DIU). The FPDE provides timing and control 
functions for the focal plane arrays located in the SFHe dewar. It consists of approximately 
six circuit boards and some associated power-regulation electronics. This unit will be 
packaged in an enclosure about 8 x 6 x 12 in. and physically located as near the arrays 
as possible. 

The DIU provides signal conditioning, control, and telemetering functions and will be installed 
in a chassis mounted within one of the sealed HH-G canisters. Enclosing the DIU in one of 
these canisters will simplify the packaging and thermal environment management tasks and 
allow use of the optical disk drive, which is only rated to 50,000-ft operating altitude. Using the 
canister will require the design and approval of a modified canister baseplate. The 
modification is necessary to add at least one additional connector to interface the DIU to the 
FPDE. Both units will be mounted with the SFHe dewar as part of a single assembly located 
on the Hitchhiker-M structure. Approximate weights of the two units are: FPDE, 40 lb; DIU, 
210 lb including the insulated HH-G canister. 

The DIU hardware will consist of three subassemblies: an erasable/rewritable optical disk 
drive, an electronics subchassis, and a power supply/filter assembly, plus some interconnect 
cable. The disk drive is a severe environment unit and will require no additional vibration or 
shock isolation in its mounting. The electronics circuitry will be arranged on seven circuit 
boards (five different types) installed in the electronics subchassis as plug-in modules. 
Circuit-board functions will be as follows: 
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(1) Central timing and sequencing controller 

(2) Disk controller and program memory 

(3) Analog-to-digital converters (three identical circuit boards) 

(4) Telemetry data encoder 

(5) Focal plane array and environmental heater control 

The DIU chassis will also contain a power supply/filter assembly, converting 28 V dc orbiter 
power to the dc voltage levels required for the disk drive and the electronics and providing 
power to the FPDE. This assembly will also contain any power devices required for driving 
environmental heaters. 

Power consumption for the total flight system will not exceed 175 W: 50 W for the FPDE and 
125 W for the DIU. There may be additional power requirements for heaters within the two 
electronics assemblies. 

The CGSE will consist of a PC AT equivalent, with keyboard, floppy and hard disk drives, 
monitor, and printer. The equipment will also include a serial data decommutator compatible 
with the medium rate downlink output from the MRGSE. The hard disk drive will be a 
SCSI-compatible high-speed device which will be capable of storing at least 300 Mbyte of 
downlink data. The CGSE will be packaged either as a portable unit or in a small, rack- 
mountable assembly. 

4.4.2.3 Operation . A block diagram of the system electronics appears as Fig. 4-5. There 
are four focal plane arrays. Each array consists of 58 x 62 elements; each array is 
completely read five times per second. Two elements in each array are sampled simultane- 
ously, producing a total of eight analog outputs at each sample time. This results in eight 
serial analog data streams clocked at an 8.99 K word/second rate. Each analog signal Is 
converted to a 12-bit digital data word and written to a digital storage medium. 

The time for one complete experiment sample, or pass, totals 15 min, with three passes 
occurring per day. The data collected form the four arrays during one pass total 64.8 M 
12-bit words, the equivalent of 97.1 Mbyte. 

The data-storage medium will be an erasable/rewritable optical storage disk. Its capacity is 
300 Mbyte, enough to store one day’s data. Since the mission time will exceed one day, 
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Fig. 4-5 Block Diagram of System Electronics 
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some means must be provided to downlink data in order to clear space on the disk. To 
perform this function, data bytes are read from the disk to the data encoder, which 
assembles them into a telemetry format compatible with the medium rate KU-band downlink; 
it also merges selected housekeeping status data and mission elapsed time (MET) with the 
array data to form the telemetry data frame shown in Fig. 4-6. 

Each telemetry frame contains 462 12-bit words (5544 bits); it takes eight telemetry frames 
to downlink the data from one focal plane array. To downlink the full amount of array data 
(four arrays) from one experiment pass will require transmitting 799.2 M bits. That is, for 
example, 20 min of downlink time at 666 K bit/s, which will be required every 8 h. 

Although the medium rate data will be collected on tape by the MRGSE, the CGSE will be 
capable of decommutating the MRGSE serial data and clock outputs in order to monitor the 
condition of the focal plane, dewar, and flight electronics systems. The CGSE will be capable 
of storing a limited amount of downlink data on disk to allow near realtime analysis of the focal 
plane data. The CGSE will also support the RS-232 synchronous up- and down-links in order 
to pass command data to the flight electronics. 

4.4.3 Preflight Tests 

The Sensor Electronics Department at LMSC's Research & Development Division has a 
complete focal-plane test and analyses facility. These facilities will be available to this 
program. Ail the requisite peripheral equipment such as test dewars and pumping stations 
are also available. A low-background test chamber with two cryogenic blackbody sources 
and two cryogenic filter wheels is under construction and will be available to the program. 
One of three FPA test stations is shown in Fig. 4-7. In this view, the FPA is cooled with an 
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open-cycle He refrigerator. A block diagram of the system is shown in Fig. 4-8. This 
general-purpose test facility has been used to characterize InSb, HgCdTe, and Si:XX FPAs. 

After delivery of the FPAs, the supplier’s FPA characteristics will be corroborated and the 
optimum operating conditions will be estimated. These tests will use the FPA test facility drive 
electronics. The response of the FPA to discrete ionizing events will be tested using an 
in-house Co 60 gamma source. Tests will be performed at dose rates low enough not to 
introduce any long-term dose effects. 

Initially, the focal plane drive electronics will be tested independently of the FPA. After 
verification of required performance levels, the electronics package will be integrated with 
the FPA and full-scale performance checks of the system will be conducted. These tests will 
include additional Co 60 gamma tests to determine the response of the system to discrete 
ionizing events. 

Similarly, the remainder of the flight electronics package will first be tested without the FPA 
and FPDE to verify performance. Finally, full system-level tests will be performed. These tests 
will include Co 60 gamma-pulse testing. 

The final tests at LMSC will be performed after the FPA system is integrated into HELD. 
These tests will include long-term thermal stability and thermal response time in addition to 
electrical performance tests. 

4.4.4 Integration With Host Spacecraft 

HELD and accompanying flight electronics will be mounted on one of the new extension 
capabilities of the Space Transport System (STS), the HITCHHIKER-M (HH-M). The HH-M, 
developed by Marshall Space Flight Center (MSFC), is a standardized mechanical platform 
which will carry up to 1 200 lb of customer equipment mounted on a cross bay bridge type 
structure in the Space Shuttle. The HH-M will also be equipped with an avionics control unit 
that gives the customer easy access to the orbiter’s resources from which a total system can 
be configured. 

HELD is shown installed in HH-M in Figs. 4-9 and 4-10. The dewar will occupy the top of two 
bays of the HH-M. The two HH-M bays can accommodate 760 lb for this configuration. 
Table 4-2 lists the HELD component weights for the system. The flight electronics box is not 
shown, but it will mount on the front of the HH-M, which can support 170 lb. 
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Rg. 4-8 High-Speed Focal-Plane Test Station 
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Fig. 4-10 HELD-HHM Integration (Side View) 

4.4.5 HELD Testing, Fluid Management, and Safety 

These subjects are described in detail in the ICE final report previously referenced. 
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Table 4-2 WEIGHT SUMMARY (lb) 


He-ll Tank and Plumbing 

106 

200 L He-ll 

55 

He-1 Tank 

9 

15 L He-1 

2 

Vapor-Cooled Shields (3) 

66 

MLI 

25 

PODS (6) 

8 

Support Ring 

52 

Vacuum Ring 

102 

External Plumbing 

5 

Miscellaneous 

10 

Total HELD Launch Weight 

440 

FPA Experiment Package 

1 

FPA Experiment Electronics 

250 

Microsphere Experiment Package 

25 

Microsphere Flight Electronics 

30 

HH-M Interface Support Structure 

50 

Total Launch Weight (lb) 

796 
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Section 5 

PRELIMINARY COST ANALYSIS 

In this section we present a preliminary cost analysis for the flight experiment Infrared 
Detector Performance in the South Atlantic Anomaly. The cost analysis is based on the 
following Work Breakdown Structure (WBS). This cost analysis addresses only those tasks 
directly related to this experiment. Testing HELD and its integration was presented in 
"Integrated Cryogenic Experiment" (NASA Contractor Report 177538). 

5.1 MANAGEMENT 

This includes the program controls on cost and schedule (including subcontractor cost and 
schedule), configuration management, administration, and technical publications. 

5.2 PRODUCT ASSURANCE 

This includes safety, reliability, quality assurance engineering, and inspection. 

5.3 ENGINEERING AND DESIGN 

This includes all the analyses necessary to complete the design of the experiment, 
specifications for focal plane arrays, design electronics, and experiment integration. 

5.4 MANUFACTURING 

This includes all the labor and materials necessary to fabricate and assemble the flight 
experiment. 

5.5 TEST 

Includes all testing, data reduction, and documentation of subsystem, system, and flight 
tests. 
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5.6 GROUND SUPPORT EQUIPMENT 

This includes design, fabrication, and assembly of all fixtures, handling equipment, service 
equipment, and containers necessary for the experiment. 

5.7 PAYLOAD INTEGRATION 

This is primarily the labor necessary to support the integration of the experiment with HH-M 
and the shuttle. 

A more detailed WBS is shown in Fig. 5-1 . The schedule for the proposed 33-month program 
is shown in Fig. 5-2. A cost breakdown by WBS task is shown in Table 5-1. The estimated 
cost of the program is $3,037,149. 
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Table 5-1 ESTIMATED MANPOWER AND COST BREAKDOWN 


1.1 

Management 

5280 hours 

1.2 

Product Assurance 

2640 

1.3 

Design 

2400 

1.4 

Manufacturing 

1440 


1.4.3 Subcontract’ 


1.5 

Test 

4640 

1.6 

GSE 

1040 

1.7 

Payload Integration 

1920 


TOTAL 

19,360 


Total Cost Excluding Fees: $3,037,149 


’For four focal plane arrays 

Significant savings accrue if detector arrays being developed for NASA programs 
are Government-furnished equipment (GFE). 
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Appendix A 

CHORD LENGTH AND PULSE HEIGHT CALCULATION 

SUMMARY 


The paper by Bradford 3 gives a general expression for the sum distribution of chord lengths, 
C (/), in an arbitrary rectangular parallelepiped of dimensions h x a x b, where h < a < b. 
The sum distribution is the probability that a chord will be of length greater than ( . The chord 
length density, c(/), the probability that a chord will have length between i and / + d/ , is 
the negative of the sum distribution with respect to the chord length, / 

°W = -^c (i) 


Bradford gives only expressions for the sum distribution, and plots the sum distribution as a 
function of chord length for various objects. 

Kellerer 2 gives general expressions for calculating both the sum distribution , C(/), and the 
probability density, c(/). These expressions are 


C(f) = 2k 



hx 2 

Vo-* 2 ) 


F(/x) + xF * (/x) 


dx 


( 1 ) 


r r 



X 3 / 


f(/x) + 2 x 2 F(/x) 


dx 


oo 



y^j 


( 2 ) 
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The second integral in c(/), expression (2), applies only when / > h, and is zero otherwise. 
Here 

d = 0 for / < h 

= /(I - h 2 // 2 ) for / > h 

F(t) = |°°f(x) dx 

F*(t) = J " F(x) dx - tF(t) 

f(t) is the two-dimensional probability density for a chord of length t in a rectangle of 
dimensions a x b. This is calculated for the specific geometry, as in Coleman, 1 and, in our 
case of a rectangular cross section, a x b, to our parallelepiped is given by 3 

f(t) = —Xu for 0 < t < a 

a + b 

_ a 2 b 

" t 2 (a + b)(t 2 - a 2 ) 1 / 2 for a < t < b 

_ ab / a b t^\ 

_ t 2 (a -I- b) y (t 2 - a 2 ) 1 / 2 (t 2 -b 2 ) 1/2 ab J 

for b < t < 7(a 2 + b 2 ) 

= 0 for 7(a 2 + b 2 ) < t 

It is from expression (1) that Bradford derives his analytical expressions for the sum 
distribution. 

Bradford gives individual solutions for the sum distribution for four regions of the chord length, 
i , related to the dimensions, h, a, and b, of the region in which the chord length distribution is 
to be found. These regions are 
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(1) 0 < / < h 

(2) h < / < a 

(3) a < / < b 

(4) b < / < /max 
where Imax = (a 2 + b 2 + c 2 ) 1/2 

There is an error in the expression for the sum distribution in region 3 in Bradford’s paper. 

Bradford's expressions for the sum distribution may then be differentiated with respect to the 
chord length to obtain analytical expressions for the probability density, c(/). This was 
attempted, but the results were disappointing due to the excessive time required to correct 
Bradford’s sum distribution expressions and to derive the probability density from them. 

Rather, we used the alternative method of calculating sum distribution and probability density 
values directly from expressions (1) and (2) by numerical integration. In this case, it is seen 
that since 


f(t) = 0 for t > y(a 2 + b 2 ) 

the upper limit of the integration in the expressions for F(t), F*(t), and c(/) becomes 

/(a 2 + b 2 ) 

The numerical integration was performed by the quartic analog to the Simpson’s Rule 
method. 

The numerical integration program gave sum distribution results which matched those of 
Bradford for the case of a cube of side 10 p. and a volume of 10 x 40 x 60 ji. The program 
was then employed to calculate the sum distribution and probability density for chord lengths 
in the 75- x 75- x 300-p. volume of interest. 

At this point, we acquired the NOVICE code, which computes pulse height distributions via a 
Monte Carlo simulation. We stopped development of the in-house chord length calculation 
model when the already maturely developed NOVICE code became available. 
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